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?H Abstract 

.^ Using the framework of transformation optics, this paper presents a detailed analysis 

of a non-singular square cloak for acoustic, out-of-plane shear elastic, and electromagnetic 

"^ waves. The generating map is examined in detail and linked to the material properties of 

I— I the cloak. Analysis of wave propagation through the cloak is presented and accompanied by 

T^ numerical illustrations. The efficacy of the regularised cloak is demonstrated and an objective 

'T^ numerical measure of the quality of the cloaking effect is provided. It is demonstrated that 

(^ the cloaking effect persists over a wide range of frequencies. As a demonstration of the 

"t^ effectiveness of the regularised cloak, a Young's double slit experiment is presented. The 

rj stability of the interference pattern is examined when a cloaked and uncloaked obstacle are 

I ^ I successively placed in front of one of the apertures. This novel link with a well-known quantum 

mechanical experiment provides an additional method through which the quality of cloaks 

may be examined. In the second half of the paper, it is shown that an approximate cloak 

may be constructed using a discrete lattice structure. The geometry and material properties 

of the lattice are derived from the continuum cloak. The efficiency of the approximate lattice 

^^ cloak is analysed and a series of illustrative simulations presented. It is demonstrated that 

T-H effective cloaking may be obtained by using a relatively simple lattice structure, particularly 

^^ in the low frequency regime. 
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%j 1 Introduction 

Following the publication of two seminal papers in 2006 by Pendry et aL pQ and Leonhardt [2], 
there has been very considerable interest and many scholarly works have been devoted to the 
concept of invisibility cloaks; for example, see the recent review article ^3j and references therein. 
The experimental validation of cloaking by Schurig et al. [1] has further increased both scholarly, 
and popular, interest in invisibility cloaks. The concept of cloaking via transformation optics is 
due to an earlier fundamental result by Greenleaf et al. [S", 1)] on singular transformations and 
applications to cloaking for conductivity. The fundamental property that enables objects to be 
cloaked from electromagnetic waves is the metric invariance of Maxwell's equations [71 [8]. 

The classical approach to cloaking via transformation geometry involves deforming a region 
such that a point is mapped to a finite region corresponding to the inner boundary of the 
cloak. Usually, such transformations involve transforming a point into a finite region with a 
smooth boundary, such as an ellipse (see [U [HI [TO! , among others) . The mapping is non-singular 



> 

in 



everywhere, except at the initial point which is deformed into the inner boundary of the cloak. 
Since the material properties of the cloak are determined by the transformation |llj . such 
singular transformations lead to material properties that are singular. Greenleaf et al. addressed 
this issue in two papers |12| IT3] and derived the condition for finite energy solutions for weak 
solutions of the Helmholtz equation and Maxwell's equations. In [12], it was demonstrated that 
finite energy solutions to the cloaking problem for the Helmholtz equation exist for an object 
with a single layer cloak. However, for the case of Maxwell's equations with internal currents, 
the cloaking of an infinite cylinder cannot be achieved with a single layer or without imposing a 
physical surface at the inner boundary of the cloak. In the same paper, Greenleaf et al. derived 
an identity linking the transformed scalar wave equation to the metric of the deformed space, 
which may then by linked to the material properties of the cloak [13] . In 2008 Norris |9j studied 
acoustic cloaking and re-derived an equivalent identity to that in |12j using the framework of 
finite elasticity, leading to a cloak with a density described by a rank two tensor. Moreover, it 
was demonstrated that the total mass of the cloak is infinite for the case of perfect cloaking. 
Norris further demonstrated that the problem of infinite mass could be overcome if both the 
density and elastic properties of the cloak were anisotropic. An alternative approach to negate 
the problem of singular material properties is to construct a so-called near cloak by regularising 
the transformation |14) . Rather than mapping a single point to the inner boundary of the cloak 
Kohn et al. proposed mapping a ball of small, but finite, radius to the inner boundary [Hj. A 
small regularisation parameter which characterises the initial radius of the ball is introduced, 
which results in a non-singular mapping on the cloak and its boundary. The regularisation 
procedure was used to create illustrative near cloaks in |9j. 

In 2006, Milton et al. [15j examined how the equations of motion for a general elastic medium 
transform under an arbitrary curvilinear transformation. It was shown that a priori requiring 
a symmetric stress tensor enforces a particular choice of the gauge (i.e. the manner in which 
the displacements transforms). It was found that, in general, the equations of motion are 
not invariant under transformation but are mapped to a more general system with non-scalar 
density. Milton et al. demonstrated that a special case of the so-called Willis equations |16j 
remain invariant under general curvilinear transformations. In [15^ identities linking the material 
properties of the transformed material to the transformation, for both classical elasticity and the 
more general Willis materials are derived. In 2011, Norris et al. [T^ further generalised the work 
of Milton et al., deriving a more general system of transformed equations without imposing the 
constraint of symmetric stress. The material properties of the transformed system were derived 
explicitly and shown to depend on both the transformation and choice of gauge. Together [15] 
and [17] provide a comprehensive framework in which to investigate cloaking in elastodynamics. 

A design for a cloak to control flexural waves in thin plates was proposed by Farhat et 
al. [H] . The cloak is constructed of several concentric layers of piecewise constant isotropic elastic 
material. Farhat et al. also presented a simplified model suitable for practical implementation 
with ten layers using six different materials. Following [TH], an experimental group led by 
Wegener fabricated a cloak based on the work of Farhat et al. using twenty concentric rings 
and sixteen different elastic metamaterials |19| . Physical measurements were compared with 
numerical simulations and found to be in good agreement. Control of in-plane waves governed by 
the Navier equations was examined by Brun et al. [20]. In [20], the authors modelled a circular 
cloak using the classical radial transformation by deforming a disc to an annulus. The efficiency 
of the cloak was illustrated using finite element simulations and the numerical solution of the 
cloaking problem was compared with the Green's function for a homogeneous elastic space. 

An influential paper f22] presents a transformation optics algorithm for a problem of electro- 
magnetism involving a cloak of a square shape. The transformation is performed in Cartesian 
coordinates and results in a piecewise smooth cloak on the interior points, with matching regions 
in the neighbourhood of corners and a singularity at the origin transformed into the inner 
boundary of the cloak. The model of such a continuum cloak received substantial attention and 
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Figure 1: (a)-(c) The field u{x) for tlie Young's double slit experiment with no inclusion, an 
uncloaked inclusion, and a cloaked inclusion respectively, (d) A plot of 1^(0;) | over the observation 
screen illustrating the interference fringes for cases (a)-(c). (An animated version of this figure may 
be found in the supplementary material [2T|.) 



subsequent use by the modelling community (see, for example, [231 [Ml US [Ml EZl ES] ) . In the 
majority of these papers, the emphasis is on the geometrical aspect of the possible shapes of 
the cloak, with examples ranging from polygonal and elliptical cloaks to heart-shaped cloaks. 
Although it is indeed interesting to see a wide range of transformations and geometries, it also 
remains important to understand the transformed problem in the context of the physical model, 
address the analysis of the transformed boundary or transmission conditions and furthermore 
derive the properties of the solutions. The paper [22], which stimulated a good level of discussion 
on the topic, also admits a deficiency regarding the analysis of the solution near the boundary 
of the cloak. Apparently, no indication is given about the sensitivity of the result to the type of 
boundary conditions (Dirichlet or Neumann) on the inner boundary of the cloak. The authors' 
evaluation of the effectiveness of the cloaking is based on a visual observation linked to a 
numerical simulation at a single frequency. On page 91, the text of [22] says "Although the 
effective fix is inaccurate in the vicinity of the boundary of the inner square, it can be shown, 
that the performance of the implemented device is not affected by this fact, which is out of the 



scope of this paper." Indeed, if the authors had attempted to change the frequency range they 
would have seen significant differences. Hence, the cloak advocated in [22] is an approximate 
cloak, where the boundary effects become important and visible as the frequency of the incident 
waves increases. 

The ideas of metric invariance in Maxwell's equations and cloaking have been taken on board 
as a technical tool and on many occasions, the researchers omit to look at the physical model 
corresponding to the transformed equations. For example, on page 99 in [26j the text reads "The 
square cloak has the same geometry as the cylindrical case, except that we replace the cylindrical 
shell by a rectangular shell with the same size" . This comparison of unlike geometries omits 
important effects, such as field concentrations near sharp corners, which make cloaking more 
difficult. Motivated by [22], Farhat et al. ^23j attempted to construct an approximate square cloak 
for out-of-plane shear waves. Farhat et al., using the method of multiple scales, introduced a 
microstructure composed of a regular array of perforations and derived a homogenised continuum 
which would approximate an ideal cloak. However, as Farhat et al. admit on page 15 of [23] "it 
is fair to say that our structured cloak is not as efficient as we would have expected" . 

Polygonal cloaks have also been the subject of experimental investigation. For example, 
Chen et al. report |29j the results of an experimental broadband hexagonal cloak based on a 
piecewise linear homogeneous transformation. Although the cloak does not render the cloaked 
object invisible, it does reduce its apparent size. The cloak is demonstrated to work for visible 
light. However, Chen et al. emphasise that the cloak only functions for light incident from 
six directions defined by the faces of the hexagon. More recently, Landy et al. produced an 
experimental uni-directional metamaterial cloak for microwaves. The reported cloak [3^ is 
based on a bilinear transformation which maps a line segment to a two dimensional region of 
space. Cloaks based on such transformations are referred to as carpet cloaks in the literature 
(see |31j among others). The advantage of such cloaks is that the requisite material properties 
are homogeneous and finite. However, Landy et al. admit that such cloaks are only effective 
over a narrow range of observation angles. The cloak is nonetheless impressive given that the 
practical implementation does not rely on the eikonal approximations as is the case with other 
implementations [32t[33l[M] . 

It appears that the work reported in [22] has generated a scope for further discussion and 
indeed further improvement of the model involving "glued transforms" that lead to approximate 
rectangular cloaks. Such cloaks are by no means exact and are frequency sensitive. A regulari- 
sation procedure, as illustrated by Kohn et al. for the spherical cloak |14j . can be applied to 
make the transformation, and hence the material properties, non-singular on the inner boundary 
of the cloak. The regularisation procedure not only simplifies the analysis, but also makes it 
physically meaningful. Furthermore, a lattice approximation is straightforward for a regularised 
square-shaped cloak. This appears to be efficient and serves a relatively wide frequency range. 

It should be emphasised that the current paper improves significantly the influential work 
by Rahm et al. [52J. In particular, the square push out transformation is regularised and it is 
demonstrated that this construction of a near cloak via regularisation does not adversely affect 
the effectiveness of the cloaking. Detailed ray analysis is carried out and the ray equations are 
derived. An objective numerical measure of the quality of the cloaking is introduced and used to 
quantify the effects of frequency, interface conditions, and orientation of the inclusion. Illustrative 
simulations are produced for cylindrical sources and the analytical ray diagrams presented. For 
physical fabrication, a discrete lattice cloak is used to approximate the regularised continuum 
transform and it is shown that the lattice structure acts as an effective cloak, particularly at low 
frequencies. 

In the present paper, the emphasis is on solutions of the Helmholtz equation. The Helmholtz 
equation arises in a wide variety of flelds including electromagnetism, elasticity, and acoustics. 
Therefore, solutions to the cloaking problem for the Helmholtz equation have a wide range of 
potential applications. However, for definiteness and ease of exposition the language of elasticity 



will be used throughout this paper. 

The current paper is structured as follows. A description of the regularised cloak in the 
continuum model of out-of-plane shear elastic waves follows the introduction. This also includes 
the discussion of the natural interface conditions on the boundaries of the cloak. An explicit 
ray tracing algorithm is developed, and the phenomenon of negative refraction on the interface 
boundaries is explained. Numerical scattering measures are included, with detailed simulations. 
The analysis also incorporates Neumann and Dirichlet boundary conditions on the inner contour 
of the cloak. As a demonstration of the effectiveness of the regularised cloak, a Young's double 
slit experiment is presented. To give the reader an idea of the results which follow, we show part 
of a simulation of the interference pattern due to a plane wave interacting with two apertures in 
figure [T] The same figure presents the fringe structure for the case of an aperture obstructed by 
an uncloaked and then cloaked inclusion. This numerical experiment convincingly demonstrates 
the high quality of the cloaking effect produced by this regularised cloak. Further discussion 
of the Young's double slit experiment is deferred to section [3} It is also shown that one of the 
undeniable advantages of such an approximate cloak is the straightforward connection with the 
discrete lattice structures. These connections are analysed in detail, and accompanied by a 
range of physical simulations. Concluding discussion and a final outline are included in the last 
section of the paper. 

2 The regularised continuum cloak 

The classical approach to cloaking via transformation geometry involves deforming a region 
such that a point is mapped to a finite region corresponding to the inner boundary of the cloak. 
Indeed, the square push out transformation proposed by Rahm et al. [22j maps a point to a 
square. The mapping is non-singular everywhere except at the inner boundary of the cloak. 
In the present paper, a regularised version of the square push out transformation is used. In 
particular, the trapezoids x are mapped to the trapezoids 0_ as illustrated in figure^ with 
continuity, but not smoothness, imposed on the interfaces between the four trapezoids. The 
mapping is non-singular on the closure of the cloak, and hence, all corresponding material 
properties are finite. It will be shown that this regularised transformation yields an effective 
broadband cloak, with finite material properties which may easily be approximated by a regular 
lattice. 

2.1 The transformation 

Consider a small square xo = {^ ■ |^i| < ^j |^2| < e} C M^, which via the map ^ is mapped 
to the square Qq = {x : \xi\ < a, \x2\ < a} C M?. Physically w is the thickness of the cloak, a 
is the semi- width of the inclusion Qq, and e is the initial semi- width of the square xo where 
< e/a <^ 1. In this case, it is convenient to decompose the cloak into four sub-domains 
X = X U . . . U X ) as illustrated in figure [2J Formally, the map ^ defines a pointwise map from 
X G x = X^^^ U . . . U x^^^ to a; = J^{X) e Q- = SlL^^ U . . . U q'lK The mapping is continuous and 
non-linear on x, and defined in a piecewise fashion such that ^ = J-^'^'{X) for X G X : where 
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aiXi + Q2 
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aiX2 + 02 
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aiX2 — 02 



with oi = w/{a + w — e) and 02 = {a + w){a — e) / (a + w — e) . The exterior of the cloak remains 
unchanged by the map, that is, X = ^{X) for X G 0+, where the bar denotes the closure of 





Figure 2: The map T maps the undeformed region x to the deformed configuration fi_. The 
boundary between £1+ and ri_ is denoted F'-*', while the interface between J7o and rj_ is denoted 
7^''. The corresponding boundaries in the undeformed configuration are denoted by F^*' and a^'^> 
respectively. 



the domain. The Jacobian matrices and determinants are then 
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It is emphasised that J^^'{xi) = detJ^^' / for a; G Q_ and e 7^ 0, that is, the map is 
continuous on both the interior and boundary of the cloak. The metric of the deformed space 



n^i^ is 5» 



(jWj»T)-i. 



The present paper will be devoted to the propagation of time harmonic out-of-plane shear 
waves of radian frequency to and displacement amplitude u(x). Lemma 2.1 in [9] allows the 
Helmholtz equation for an isotropic homogeneous medium fiVx • C^ x)u{X) + guj'^u{X) = for 
X E X to be written in deformed co-ordinates as 



[V • (C(*) (a;) V) + p(^) {x)uj'^]u{x) = 0, a; e fi 



(i) 



(1) 



where fi is the constant ambient stiffness, g is the constant ambient density, C^'^'{x) = 
[n/ J^'^' {x)]J^^\x)[J^^' (x)]"^ is the transformed stiffness matrix and p^^'{x) = q/J^\x) is the 
transformed density. The differential operator Vx is written in the undeformed space and 
should be distinguished from V which is written in the deformed coordinates. 

Since the mapping is continuous on f2_, the material properties of the cloak are non-singular. 
The transformed stiffness tensor is symmetric and positive definite. The density is a scalar and 
bounded on VL-. Physically, the transformed material properties correspond to a heterogeneous 
anisotropic medium. 

2.2 Interface conditions 

Without loss of generality, it is convenient to restrict the following analysis to a single side of 
the cloak. With reference to figure 2l consider the sub-domain 0_ C M^ in the absence of the 
inclusion and remaining three sides of the cloak. In the absence of sources the out-of-plane shear 
deformation amplitude of an outgoing time-harmonic wave of angular frequency u satisfies the 
following equation 

Cu{x) = 0, (2) 

together with the Sommerfield radiation condition at infinity. Here £ = V • (A(ic)V) + p{x)uP 
is the Helmholtz operator, A{x) and p{x) are defined as 

I pi for X € r2+ I Q for x € 1^+ 

Let v{x) be a continuous piecewise smooth solution of the Helmholtz equation in M^ satisfying 
the Sommerfield radiation condition at infinity. Integrating the difference u{x)Cv{x) — v{x)Cu{x) 
over a disc Dr of radius r containing il_ yields 



/<■ 



= / fnV • AVv - vV ■ AVu) da;, 

(u~n ■ AVv~ — v~n ■ AVu~) dx — I (u'^n ■ AVv"*" — v^n ■ AVu'^) dx 



an*'' an''' 



+ p {un ■ Vv + vn ■ Vn) da;, 

where the fact that Vn • AVv = Vv ■ AVu (since A is symmetric) has already been used. Since 
u{x) and v(x) represent outgoing solutions, the final integral vanishes as r — ?■ oo. Thus, the 
essential interface condition is continuity of the field 

[u]=0 on dn^y, (4) 

and the natural interface condition is continuity of tractions 

n-C'^'^^Vu- = pn-Vu+ on dn^ . (5) 

These interface conditions are equivalent to those derived in [9J for acoustic pressure by application 
of Nanson's formula and imposing particle continuity. 



2.3 The cloaking problem 

Consider the propagation of time harmonic out-of-plane deformations, generated by a point 
source, in a homogeneous infinite elastic sohd in which is embedded an inclusion surrounded by 
a cloak. The displacement amplitude field then satisfies 

[V ■{A{x)\/)+ p{x)uj'^]u{x) = -5{x-xo), xgR^X^o, xq £ n+ (6) 

[/ioV • (V) + Qouj^]u{x) = 0, x€ Qo, (7) 

with continuity of u{x) and tractions on all internal boundaries according to Q and ([5|. 
Additionally, the Sommerfeld radiation condition is imposed at infinity. The stiffness tensor 
A{x) and density p{x) are 

I fil for X G il+ 1 ^ for £C G r2_|_ 

and fiQ and po are the stiffness and density of the inclusion respectively. 

2.4 The ray equations 

Whilst, in principle, the exact wave behaviour of the displacement field can be obtained from 
the solution of the cloaking problem, it is illuminating to consider the leading order behaviour of 
rays through the cloak. Consider a WKB expansion of the displacement amplitude field in terms 
of angular frequency w, and the amplitude and phase functions Un{x) and ^p{x) respectively 



oo 



u{x) ^ e'"'^(") Yl '"^"i"^^ ' as^^oo, (9) 

n=0 

whence the leading order equation for the phase on the interior of the cloak is 

H{x,s) = 0, (10) 

where H{x,s) = pg^^s ■ g ^s — 1, s = \7(p is the slowness vector, p and g are the stiffness 
and density of the ambient medium respectively, and g is the metric of the transformation. In 
terms of wave propagation, the conserved quantity H{x, s) represents the first order slowness 
contours |35j . The characteristics of the quantity H{x, s) then satisfy the following system 

dH _ dx _dH ds _ dH 

where t is the ray parameter (equivalently, the time-like parameter). At this point, it is convenient 
to introduce index summation notation where summation, from 1 to 2, over repeated indices is 



implied. The system (|11|) may then be expressed as 

-2g^^pSmSnJnl-^, ^ = '^Q^^ fJ-JilJjlSj , (12) 



dt oxj dt 



where Jij = {J)ij are the components of the Jacobian matrix and should be distinguished from 
the J, the Jacobian determinant. The superscript labels have been omitted for brevity, but Jij 



j{k) 

the cloak. 



and J should be understood as J^ • and J^^> for /c = 1, ... 4 corresponding to the four sides of 



From equation (10), an alternative representation, in terms of wave normals n and the phase 
velocity v, is 

pg^^n ■ g^^n — u^ = 0. (13) 



The representation (13) is obtained by assuming a plane wave solution to the Helmholtz equation 

Alternatively, seeking a solution of the full wave equation in the form of 



(see, for example [35] 

the leading term in a WKB expansion 

oo 
n=0 






as uj 



oo, 



yields the same result with di^/dt = v. Prom (10) and (13) the slowness vector can be expressed 



in terms of the original material properties (through g and /i) and the map (through J) as 



n 




(14) 



Further, in the undeformed configuration, the equivalent conserved quantities are jj-q S-S—1 = 



and fig ^ =V'^. Together with (10) and (13), these two equations imply that 



J'S 



J-^N 



V 



J-'^N 



(15) 



Now, consider a ray (line) in the ambient medium, in direction N passing through Xq and 
parameterised by t. The corresponding curve in the cloak is x(t) = ^{Xq + tN), whence 







dxi 
dt 


= J^JN,. 


which 


using (15) may 


be rewritten thus 








dxi 

dt 


JilJjlSjJ 


Taking 


the derivative 


of (15) for constant N yields 



(16) 



dsj 
dt 



SkS-n'Jkj'Jlm'-'n 



dJ- 



j« 



dxi 



Using the compatibility condition that the deformation gradient should be irrotational under 

finite deformation ejkdJ~^ik/dxj = Oj, the partial derivative above may be written as dJ^^ji/dxi, 

whence 

dsi _ _ dJmi rjj^ ..„^ 

dt OXi V Q 



where Sjk is the permutation tensor and the equality JimdJji /dxi = — J-^ dJim/dxi has already 
been used. Consider the characteristic equations for the waves in the cloak 



the equations of the transformed rays (16) and (17). The system (16) and 




together with 
are the equations 

of characteristics in the cloak, up to an arbitrary scaling constant of 2-v//u/^. Thus, to leading 
order, rays (or straight lines) in the ambient medium map directly to rays in the cloak. 

Figure [3] shows rays emanating from a point source, passing through the cloak and emerging 
from the cloak on their original trajectory. In this sense, the object is "invisible" to an observer 
outside the cloaking region. The figure clearly illustrates how wave propagation in the cloak is 
related to the map. Animated versions of figureplcan be found in the supplementary material |21| . 

An interesting alternative perspective is apparent if figure [3] is viewed, not as rays diverging 
from a source, but as rays converging to a focal point. It is observed that the rays converge to 
the focal point around the inclusion. One can envisage several applications where such an effect 
may be useful. For example, image distortion from the mirror mounts in telescopes could be 
reduced by cloaking the mounts. In addition, apparatus and mounting structures on microwave 
receivers could be cloaked to improve the quality of the signal. One could also conceive of 
cloaking mounting points and the surrounding structures in laser cutting machines to protect 
them from accidental damage. 
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Figure 3: Plots of the ray paths through the cloak for a cylindrical source. The grey lines indicate 
the deformation of the space inside the cloak. (Animated versions of these figures may be found in 
the supplementary material |21j.) 



2.4.1 Negative refraction 

It is evident from figure [3] that whilst the rays are continuous, they are not necessarily difFer- 
entiable. In particular, at the interface between the cloak and ambient medium, and at the 
internal interfaces of the cloak, refraction occurs characterised by a discontinuity of the spatial 
derivatives of the rays. Of particular interest are the regions on the outer boundary of the cloak 
where negative refraction occurs. 



Consider figure 3a Negative refraction occurs on the right hand interface between the cloak 
and the ambient medium. A ray exiting the right hand side of the cloak with gradient M at 
point X(°) = a;(0) may be described by the equation X^^'^ - xf ^ = M{x['^ - xj°^) in the 
ambient medium, where X^^' is the position of the source. The behaviour of the ray at the 
interface is entirely characterised by the position of the source relative to the interface and its 
initial gradient. Therefore, the following analysis may be restricted to the right hand side of the 
cloak. On the interior of the right hand side of the cloak, the ray exiting the cloak at a^''^^ is 
characterised by 



.(0) 



{0)^ 



X2 



xi\ M + 



Xl — Q2 



The gradient of the ray as it approaches the exterior boundary from the interior of the cloak is 

then 

^ , dx2 M(a + w — e)(a + w) — x), '(a + e) 

m = hm - — = — ^ —. ^r ^-^ -. 

xi^{a+w)- dxi w[a + w) 

Thus, the gradient is discontinuous at the exterior interface. For negative refraction it is required 
that m*M < 0, which leads to the following inequalities 



< M < X, 



(0) 



a + e 



(a + w){a + w — e) 



or X. 



(0) 



a + e 



{a + w){a + w — e] 



< M < 0. (18) 



Assuming the source lies on the line X2 = as in figure [3a| the inequalities reduce to the single 
inequality 



yis) ^ {a + w){w-2e) 
^ a + e 



10 



which is satisfied for all sources outside the cloak X}* < —{a + w), if w < a + 3e. Thus, for a 
sufficiently thin cloak and a cylindrical source placed along X2 = and any distance from the 
cloak, negative refraction is expected on the opposite side of cloak. 



Alternatively, for a source located along the line Xi = the inequalities (18) become 



^ (a + w-e) {a + w- e) ^ ^ ^ 

where the fact that \x2 \ < (a + w) has been used. Since a,w > 0, and < e ^ 1, the 
above inequalities are never satisfied. Hence, the lack of negative refraction on the horizontal 



interfaces in figure 3a Similar arguments may be used to consider other regions where negative 
refraction may, or may not occur. It is observed that negative refraction always occurs at the 
interface between the four regions of the cloak, where the material properties (equivalently the 
transformation) are not smooth. 

2.4.2 Scattering measure 

It is desirable to have some quantifiable measure of the quality of the cloak with respect to 
shielding, rather than relying on visual observations. However, it is not obvious what "quality" 
means with respect to a cloak, given that there are essentially three fields involved, i.e. the ideal 
field in the absence of both cloak and inclusion, the uncloaked field with an inclusion present but 
without a cloak, and the cloaked field with both the inclusion and cloak. Previous experimental 
works [19] have used an L2 norm computed directly from the measured fields to place a numerical 
value on the quality of the cloak. It is in this spirit that the following ^^ scattering measure'^ is 
formally introduced as a tool to quantify the cloaking effect 



£{ui,U2,n)=\ \ui{x) - U2{x)\^ dx \ \u2{x)\'^dx\ , (20) 




where 7?. C M^ is some region outside the cloak, and ui[x) and U2{x) are any two fields. In 
the present paper, the quantities £{uu,uo,TZ) and £{uc,uo,TZ) are presented for a series of 
illustrative simulations. The field uq{x) = ^Hq {oJ^J~pJ]x\x — xq\)/4, is the Green's function for 
the unperturbed problem and represents the "ideal" field, Uu{x) and Uc{x) are the uncloaked and 
cloaked fields respectively. Thus, perfect cloaking corresponds to a vanishing £. Along with the 
raw scattering measures an additional quantity, Q = \£{uu-,uq,'R,) — £{uc,uo,TZ)\/£{uu,uo,TZ) 
characterising the relative reduction of the scattering measure by the introduction of a cloak is 
also presented. It should be emphasised that this is only one of a number of possible measures 
of quality. 

Choice of TZ. For the purpose of illustration three different regions of integration are considered 
as shown in figure Hi The three regions used were chosen as follows: (a) TZi is the most strict 
region used taking into account significant near field effects and a wide range of scattering 
angles. However, it is unlikely that this region would be measurable in practice, (b) The forward 
scattering region (7^2) is relevant if the scattered field is measurable over a wide range of forward 
scattering angles, (c) The corner scattering region (T^s) is employed for sources located along 
the diagonal of the square inclusion. It is emphasised that ||7^i|| 7^ 1 1 '^2 1 1 = ll'^sll- 

In the following section the scattering measures will be presented for a series of illustrative 
simulations. 

2.5 Illustrative simulations 

A series of illustrative simulations were created using the finite element software COMSOL 
Multiphysics® . For the purposes of these computations, the following non-dimensional parameter 
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(a) Strict scattering region (b) Forward Scattering (c) Corner Scattering 

Figure 4: The three regions used for computation of the scattering measure. 



Source 
Position Frequency 



Scattering Measure £ 
Uncloaked Cloaked Q 



[-3,0]T 
[-3,0]T 
-3,3]T/V2 
-3,3]T/y2 



Scattering region IZi 
5 0.1529 

10 0.1455 

5 0.2002 

10 0.3286 



4.351 X lO"'' 
4.514 X 10-"^ 
3.941 X 10""^ 
4.068 X 10-"^ 



0.9972 
0.9969 
0.9980 
0.9988 







Scattering region 7^2 




[-3,0]T 


5 


0.3224 3.664 x 10"^ 


0.9989 


[-3,0]T 


10 


0.3093 1.167 X 10-3 


0.9962 






Scattering region TZ^, 




[-3,3]T/V2 


5 


0.2988 3.654 x 10"^ 


0.9988 


[-3,3]T/V2 


10 


0.2988 7.803 x 10"^ 


0.9974 



Table 1: The scattering measures corresponding to the simulations shown in figures p^ and pi 

values J were chosen: a = 0.5, w = 0.5, fi = g = 1, fio = 0.1, ^o = 0,e = 1 x 10~^. Figures [5 
and [6 show the displacement field u{x) for a cylindrical source oscillating at a; = 5 and a; = 10 
respectively. The figures clearly illustrate the efficacy of the square cloak, even at relatively high 
frequencies. Table [T] shows the corresponding scattering measures as introduced in section 2.4.2 



It is clear that this square ^^push out^^ cloak is highly effective. Indeed, for the illustrative 
simulations presented here, the cloak reduces the scattering measure by not less than 99.62% 
compared with the uncloaked inclusion. 

Figure [7] shows the scattering measure plotted against non-dimensional angular frequency uj 
(with ^ = ^ = 1). The solid curve in figure [7] corresponds to the continuum, in the absence of 
both cloak and inclusion. This curve gives an indication of the numerical error in the simulation 
induced by, for example, the use of perfectly matched layers and the numerical discretisation. 
The dashed curve corresponds to the cloaked inclusion, whilst the dotted curve corresponds to the 
uncloaked inclusion. It is observed that the numerical measure of the cloaked inclusion remains 
close to that of the intact continuum for a large range of frequencies. Moving to dimensional 
quantities, suppose the simulation corresponded to a particular polarization of an electric wave 



^Throughout this paper, all numerical parameters are normalised such fj, = g = 1 unless otherwise stated. 
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(d) Cloaked, xq = [-3, 3]'^/V2 

Figure 5: Plots of the field u for the uncloaked and cloaked square inclusion, where the angular 
frequency of excitation is w = 5. The position of the source is indicated under the relevant plot. The 
colour scale is linear from blue (minimum) through white (zero) to red (maximum). 



travelling through glass at a speed of approximately 2 x 10^ m/s. The line cu 
then corresponds to a frequency of approximately 340 MHz. 



10 on figure [t] 



2.5.1 Boundary considerations 

Whilst cloaking via transformation geometry has been extensively treated in the literature, the 
sensitivity of the cloaking effect to the boundary conditions is rarely discussed. The cloak is 
formed by deforming a small region (a point in the case of the classical radial transformation [Ij), 
into a larger finite region. If the region is an inclusion, then the natural interface conditions may 



be determined following the method outlined in section 2.2 If the cloaked region is a void or 
rigid inclusion, however, there is some freedom in choosing the boundary condition, subject to 
the constraints of the physical problem. Figure p] shows the field u{x) for a cloaked void, with 
Neumann (parts (a) and (b)) and Dirichlet (parts (c) (d)) conditions applied to the interior of 
the cloaked region. The corresponding scattering measures are shown in table [2] 

Although the square cloak is effective in both cases, it is clear from both the figures and 
the table of scattering measures that the type of boundary condition imposed on the cloaked 
object affects the quality of the cloaking. Indeed, for a void (Neumann) the cloaking reduces 
the scattering measure by between 99.7% and 99.9% for both cj = 5 and u = 10. In contrast, 
cloaking reduces the scattering measure of a rigid inclusion (Dirichlet) by between 96.5% and 
96.8% for w = 5 and between 96.9% and 96.7% for oj = 10. The effect of the boundary condition 
may be interpreted in the following way. As a result of the transformation, the cloaked object 
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(c) Uncloaked, a;o = [-3, 3]T/y2 (d) Cloaked, ccq = [-3, 3]'^/y2 



Figure 6: Plots of the field u for the uncloaked and cloaked square inclusion where the angular 
frequency of excitation is w = 10. The position of the source is indicated under the relevant plot and 
the inclusion is located at the centre of the image in all cases. The colour scale is linear from blue 
(minimum) through white (zero) to red (maximum). 

and cloak together behave as if the void is small. In this sense, the cloaked inclusion represents 
a singular perturbation of the fundamental solution of the Helmholtz equation. In the case of a 
free void with Neumann conditions, the leading order term in the asymptotic expansion is the 
dipole term, which is of order e^ and decays like the first derivative of the fundamental solution. 
On the other hand, for a fixed void with Dirichlet conditions, the leading order term in the 
expansion is the monopole term which is of order e and decays like the fundamental solution. 
Thus, the perturbation from the free void is smaller than the perturbation from the fixed void, 
leading to improved cloaking. 

3 Cloaking path information 

In recent years there has been much interest in experiments to elucidate the fundamental 
principles of quantum mechanics, and in particular the relationship between measurement 
and system behaviour. One basic experiment which with its variants features in many such 
experimental studies is the classical Young's double slit experiment (see, for example, [36]). This 
suggested to us that it would be of interest to consider the interaction of the excellent mechanical 
cloaking demonstrated earlier with the foundational quantum mechanics experiment. 

Thus, a Young's double slit experiment is considered where a monochromatic plane wave is 
incident on a screen with two apertures. Due to the superposition of the waves passing through 
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Figure 7: (a) The figure shows the the scattering measure plotted against angular frequency, (b) 
The figure shows the log of the scattering measure plotted against angular frequency. The solid 
line corresponds to the continuum in the absence of both an inclusion and cloak. The dashed line 
represents the cloaked inclusion and the dash-dot line corresponds to the uncloaked inclusion. The 
region TZy (see figure Inland the associated text) was used to compute the error measure. 



Source 




Scattering Measure E 




Boundary Condition 


Frequency 


Uncloaked 


Cloaked 


Q 




Scattering 


region Tli 






Neumann 


5 


0.1624 


4.351 X 10""^ 


0.9973 


Neumann 


10 


0.1558 


4.540 X 10-"^ 


0.9971 


Dirichlet 


5 


0.2931 


1.038 X 10-2 


0.9646 


Dirichlet 


10 


0.2553 


7.875 X 10-3 


0.9692 




Scattering 


region 7^2 






Neumann 


5 


0.3436 


3.664 X 10-^ 


0.9989 


Neumann 


10 


0.3258 


1.163 X 10-3 


0.9964 


Dirichlet 


5 


0.4864 


1.566 X 10-2 


0.9678 


Dirichlet 


10 


0.5030 


1.673 X 10-2 


0.9667 



Table 2: The scattering measures for a void with Neumann and Dirichlet boundary conditions. 
Here the source is located at [—3,0]"^. 



the two apertures, the distinctive double slit interference pattern is produced on an observation 
screen placed on the opposite side of the apertures to the source. The result of a simulation 



of the standard experiment is shown in figure la, with the diffraction pattern produced on 



the observation screen (in this case, a vertical line near the right hand edge of figures la 



shown as curve (a) in figure Id, Placing an object (inclusion) over one slit, as in figure 



Ic) 



lb, 



partially destroys the diffraction pattern. The corresponding pattern on the observation screen 



is shown as line (b) in figure Id However, coating the object with the square 'push out cloak 
presented earlier, as shown in figure [Tel restores the original diffraction pattern almost entirely. 



The interference pattern corresponding to the cloaked object is shown as curve (c) in figure Id 



The simulation, shown in the supplementary material, confirms that the excellent cloaking 



for the inclusion position of figure [Tc| exemplified in figure Id, holds irrespective of the inclusion 
position. It has thus been conclusively demonstrated that the cloaking is of sufficient quality 
to render the interference pattern almost immune to movement in the position of the cloaked 
obstacle. In particular movement of the cloaked obstacle, it would seem, does not yield any 
information about the passage of waves through one slit or the other. This consideration would 
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(a) Uncloaked, Neumann 





(c) Uncloaked, Dirichlet 



(d) Cloaked, Dirichlet 



Figure 8: Plots of the field u for the uncloaked and cloaked square inclusion with Neumann 
boundary conditions in parts (a) and (b), and Dirichlet boundary conditions in parts (c) and (d). 
Here the source is located at a; = [—3,0]'^ and oscillates at a; = 10. The colour scale is linear from 
blue (minimum) through white (zero) to red (maximum) 

be important if one were able to carry out an experiment in which single quantised elements of 
flexural vibration were in the system at any given instance in time. The quantum mechanical 
view would be that, if no path information were available from measurements, the interference 
fringes behind the double slit should persist. 

This proposed quantum experiment raises interesting questions if an appropriate vibration 
transducer were embedded within the cloak, so that information about vibrations moulded by 
the cloak were available to experimentalists. One would assume, in line with the results of say 
optical experiments of the type referred to in [36] , that any path information gained in this way 
would be evident in a change in the fringe pattern. This suggests the interest of a comprehensive 
quantum mechanical treatment of the interaction between mechanical cloaks and measurement 
systems. 

4 Cloaking with a lattice 

Cloaks designed using transformation optics may have such extreme physical attributes that 
the requisite materials cannot be physically realised without recourse to metamaterials. It is 
with this motivation in mind that the following approximate cloak in the low frequency regime 
is developed. The cloak is constructed as an approximation to the continuum square cloak 
considered earlier, but is realised using a discrete lattice structure, formed from rods and point 
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masses. The advantage of a discrete structure over a continuous material is that much higher 
contrasts in material properties are easily realisable using lattices. The development of an 
approximate cloaking material using a lattice may allow the practical construction of cloaks. In 
the following discussion, it is emphasised that repeated indices are not summed over. 




Figure 9: 

cloak. 



The lattice formed from the principal directions of the stiffness matrix for the continuum 



With reference to the formulae for the Jacaobian of the transformation in section 2.1, the 
symmetric stiffness matrices C^"^' = [/u/ j'*^] j'*^ J^'^' are positive definite. Therefore, the stiffness 
matrix admits the following diagonalisation 



(j{i) — p{i) \^{i) p(i) 



(21) 



where P^^ = [e(\ef] 

(eigenvectors) of C'*\ and A^*^ = diag(A^ , Ag ) is the diagonal matrix of the corresponding 



are the matrices with columns consisting of the principal directions 



ordered [positive] eigenvalues such that A^ > A, 



i« 



The eigenvectors yield the principal lattice 

In 



vectors of the locally orthogonal lattice with homogenised stiffnesses A^- in direction e- . 
particular, the lattice nodes lie at the intersection points of the solutions of the following 
non-linear system of first order differential equations 



d (j) 
dr ■' 



.(■Ofji) 



X- 



+ 4". 



for i = 1, ... 4, and j = 1, 2, 



(22) 



where x- is the position vector along the characteristic defined by e- inside the i side of the 

cloak, T parametrises the curve, and x j is the initial position. Naturally, this would lead to a 
lattice with curved links. However, for a sufficiently refined lattice the curved members may 
be replaced with linear links. The lattice links are then the linearisation of the characteristic 
between two neighbouring nodes on the characteristic. Figure [9] shows the geometry of the 
lattice formed from the principal vectors of the stiffness matrix. Requiring local conservation 
of flux allows the stiffness of the lattice link parallel to e- to be determined as iijX- , where 

(i) 

iij is the length of the link along e^ . The distribution of nodal mass may be determined by 
evaluating the integral 



m{xp) 



p{x) dx, 



A{xp) 



over the unit cell A{xp) containing the lattice node at Xp. 

In principle, the lattice cloak may be constructed exactly as described above and illustrated 
in figure^ However, for narrow cloaks where w/a <^ 1, the locally orthogonal lattice depicted 
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in figure |9] may be approximated by a globally orthogonal regular square lattice. A regular 
square lattice is more convenient to implement compared with the non globally orthogonal lattice 
generated from the eigendecomposition of the stiffness matrix. Although the geometry of the 
approximate lattice is regular, it should be emphasised that the stiffness of the links and mass of 
the nodes vary with position according to the projection of A{x) and p{x) as described above. 

4.1 Geometry and governing equations for an inclusion cloaked by a lattice 

Consider a square inclusion Qq = {x : \xi\ < a, \x2\ < a}, a > 0, embedded in M^, surrounded 
by a cloak ft- = {x : a < \xi\ < a + w,a < \x2\ < a + w} \ilo, where ti; > is the thickness 
of the cloak. The cloak consists of a discrete lattice structure with lattice points at a; = ip, 
where p £ I? H {n : £n G Q^}. The lattice is statically anisotropic with links parallel and 



perpendicular to the boundaries having contrasting material properties, as shown in figure 10 




Figure 10: The lattice cloak $!_, surrounding the square inclusion JIqi embedded in the ambient 
medium 0+ . The tliick black lines in the lattice cloak indicate links of high stiffness or conductivity, 
while the thick grey lines indicate links of low or stiffness conductivity. 

As for the continuum cloak, solutions of the Helmholtz equation are of primary interest. In 
particular, the following problem for the field u{x) is studied 

[^V • (V) + guj'^]u{x) = -5{x 

[/ioV • (V) + Qouj^]u{x) 
m{p)uj^u{p) + Y^ (■ri{q,p)[u{p + q) 

where e^ = [dii,6i2\^, p ^T? , and J\f = {±ei, ±62} is the set of nearest neighbours. The stiffness 
of the lattice links are the restriction of the eigenvalues of the stiffness matrix to the links. In 
particular, for the link connecting nodes p and p + q, ri{q,p) takes the value ^i \[£p/(p+q)] if the 
vector q is parallel to the exterior boundary of the cloak, r^*^ and ^2 \[ep/{p+q)] otherwise. The 

Here, ^j \e(q-p) indicates the restriction 



Xq), X,Xo € Q.+ , 


(23) 


0, X£ fto, 


(24) 


j)-u{p)] = 0, inO_, 


(25) 



corner regions are matched as illustrated in figure 
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of A^ to the line [Ip, £{p + q)] . The associated interface conditions corresponding to continuity 



of tractions are 



n • Vu{x) 



for a; G dO, and x ± iq ^ ^_ 

r]{^q,p)u{x±eq)/n for a; € Ui T^ and £C ± q € f7„ , i = 1, . . . ,4, (26) 

^r]{^q,p)u{x ± iq)/lJ-o for x € Uj 7*"*^ ^ind x± q e $!_ 



and the Sommerfeld radiation condition at infinity. The quantity ri(q,p) is the projection of the 
diagonahsed stiffness matrix onto the lattice link connecting lattice points p and p + q. 



Physically, (23)-(26) corresponds to the problem of the propagation of time-harmonic waves 
of angular frequency u generated by a point load at xq. The field u{x) then corresponds to the 
out-of-plane displacement amplitude and g is the scalar density. The region J7_ consists of an 
array of nodes of mass m, connected by massless rods of length i and stiffness according to their 
orientation. 

4.2 Illustrative lattice simulations 

The approximate lattice cloaks were examined using the finite element software Comsol Multiphysics^ 
PML conditions were imposed on the boundary of the computational domain in order to simulate 
an infinite domain. For the purpose of illustration, a square of semi-width a = 0.5, surrounded 
by a lattice cloak with w = 0.1 and links of length 5 x 10~^ was used. The inclusion is located 
at the origin of the computational window. 

4.2.1 A basic lattice cloak 

Before proceeding to the illustrative simulations for the regular lattice with heterogeneous 
distributions of stiffness and mass, it is instructive to consider a simple approximation. Many 
cloaks created from transformation optics have the general characteristic of having a high phase 
speed parallel to the boundary of the cloak, and a phase speed in the direction normal to the 
boundary (see [37] among others). Therefore, as an initial approximation, the case of a regular 
square lattice with a homogeneous, but orthotropic distribution of stiffness and a homogeneous 
distribution of mass is considered. Consider the right-hand side of the cloak $!_ . For a narrow 
cloak with w/a <^ 1, xi ^ a-\- w and hence the density may be approximated by /? ~ 1 + a/w. 
The greatest contrast in stiffness occurs at X2 = 0, thus the vertical links are assigned stiffness 
Xl {a+w,0) and the horizontal links stiffness Ag {a + w,0). The mass of the nodes is i'^{l + a/w). 
The material properties of the remaining three sides of the cloak are adjusted accordingly. 



Figures 11 and 12 show the field u{x) for the uncloaked inclusion (a) and (d), and the 



inclusion cloaked with this basic cloak (b) and (e) . For u; = 3 figure 1 1 indicates that the basic 



cloak partially mitigates the shadow cast by the inclusion and acts to reform the cylindrical wave 



fronts behind the inclusion. As illustrated by figure 12, this partial cloaking effect deteriorates 
with increasing frequency. Indeed, in some cases, the presence of the lattice cloak seems to 
increase the shadow region. Table [3] details the values of the scattering measures for the fields 



illustrated in figure 11 and 12, The scattering measures shown in table [3] suggest that, although 



visually the basic lattice cloak appears to work reasonably well, this may not be the case. The 
fact that the basic lattice cloak increases the scattering measure compared with the uncloaked 
inclusion further emphasises the need for an objective measure of the quality of cloaks, rather 
than simply relying on visual observations. 

This increase in the scattering measure by the basic lattice cloak motivates the introduction 
of the following refined model. 
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Source Scattering Measure £ 

Position Frequency Uncloaked Cloaked Q 

Scattering region TZi 



[-3,0]T 


3 0.1430 


0.1662 


0.1617 


[-3,3]T/V2 


3 0.1113 


0.1816 


0.6327 


[-3,0]T 


5 0.1529 


0.2495 


0.6318 


[-3,3]T/V2 


5 0.2002 


0.3538 


0.7676 




Scattering region 7?-2 






[-3,0]T 


3 0.2341 


0.3362 


0.4363 


[-3,0]T 


5 0.3224 


0.4671 


0.4489 




Scattering region TZ^ 






[-3,3]T/V2 


3 0.1578 


0.3455 


1.189 


[-3,3]T/V2 


5 0.2988 


0.6011 


1.012 



Table 3: The scattering measures corresponding to the simulations for the basic lattice model 



shown in figures 11 and 12 



4.2.2 A refined lattice cloalc 



Consider now the lattice described in section 4.1 i.e. the regular square lattice with inhomoge- 



neons distribution of stiffness and mass. Figures 11 and 12 show the field u{x) for the uncloaked 
inclusion and the inclusion with a lattice cloaking. With reference to the simulations for the 
basic cloak (b) and (e) the refined lattice cloak (c) and (f), it is observed that the efficiency of 
the refined lattice cloak, whilst not as high as that of the continuum cloak, is much greater than 
that of the basic cloak. The table of scattering measures for the approximate cloak is shown in 
table |4] and further evidences the effectiveness of the refined lattice cloak. Indeed, for several 
simulations (in particular those where the scattering measure is taken over the forward or corner 
scattering regions IZi and IZ2 respectively) the efficiency of the refined cloak in reducing the 
scattering measure approaches that of the continuum cloak. 

As expected the effectiveness of the lattice cloaks reduce with increasing frequency. However, 
for sufficiently low frequencies the refined lattice cloak in particular, works well. 

5 Concluding remarks 

This work represents a comprehensive treatment of a non-singular cloak for a square inclusion. 
The significant advantage of this continuous cloak is the straightforward correspondence with a 
discrete metamaterial lattice structure. Such a connection may present a method through which 
a physical cloak may be fabricated. The material and geometric properties of the discrete cloak 
are directly linked to the properties of the continuum cloak, and hence, to the properties of 
the formal map. The effectiveness of such discrete cloaks, particularly at low frequencies, was 
demonstrated through numerical simulations and the use of objective scattering measures. 

Particular attention was paid to the objective measurement of the quality of the cloaking 
effect. The quality of the cloaks were primarily assessed using a scattering measure introduced 
as an L2 norm of the difference between the cloaked field and the ideal unperturbed field. A 
further demonstration of the efficacy of the square push out cloak was presented via the classical 
Young's double slit experiment. It was shown that the interference pattern on the observation 
screen was heavily modified when an obstacle was place in front of one of the apertures. However, 
if the obstacle was cloaked then the interference pattern remained almost entirely unperturbed. 
This numerical experiment presents a further, perhaps more interesting, method through which 
the quality of particular cloaks may be examined. Moreover, the experiment raises interesting 
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(a) Uncloaked 



(b) Basic cloak 



(c) Refined cloak 




(d) Uncloaked 



(e) Basic cloak 



(f ) Refined cloak 



Figure 11: Plots of the field u{x) for a cylindrical wave incident on a square inclusion in the 
absence of a cloak (parts (a) and (d)), a square inclusion coated with the basic lattice (parts (b) and 
(e)), and an inclusion coating with the refined lattice (parts (c) and (f)). Here the angular frequency 
of excitation is w = 3 and the source is located at Xq = [—3, 0]'^ in (a)-(c), and at Xq ~ [—3, 3]'^/\/2 
in (e)-(f). The colour scale is linear from blue (minimum) through white (zero) to red (maximum). 

questions regarding the interaction between cloaking and quantum mechanics. 
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(a) Uncloaked 



(b) Basic cloak 



(c) Refined cloak 




(d) Uncloaked 



(e) Basic cloak 



(f ) Refined cloak 



Figure 12: Plots of the field u{x) for a cylindrical wave incident on a square inclusion in the 
absence of a cloak (parts (a) and (d)), a square inclusion coated with the basic lattice model (parts 
(b) and (e)), and an inclusion coating with the refined lattice (parts (c) and (f)). Here the angular 
frequency of excitation is w = 10 and the source is located at Xq = [—3,0]'^ in (a)-(c), and at 
Xq = [—3,S\'^/\/2 in (e)-(f). The colour scale is linear from blue (minimum) through white (zero) to 
red (maximum). 
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